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Abstract 

We consider the asymptotic behaviour of the solution of one di- 
mensional stochastic differential equations and Lagevin equations in 
periodic backgrounds with zero average. We prove that in several 
such models, there is genrically a non vanishing asymptotic velocity, 
despite of the fact that the average of the background is zero. 



1 Introduction. 

We consider the one dimensional diffusion problem 



dtu = ( ^d^u + b{t, x)u ) (1^ 
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where 6 is a regular function periodic of period T in the time variable t G 
and of period L in the space variable x G R. This equation and related 
equations discussed below appear in many natural questions like molecular 
motors, population dynamics, pulsed dielectrophoresis, etc. See for example 
[6], [5], [To] [13], [E] and references therein. 

We assume that the initial condition u{0, x) is non negative and of integral 
one with |x|m(0,x) integrable, and denote by u{x,t) the solution at time t. 
Note that the integral of u with respect to x is invariant by time evolution, 
the integral of xu{x, t) is finite. 

One of the striking phenomenon is that even if the drift has zero average, 
this system may show a non zero average speed. There are many results on 
homogenization theory which can be applied to equation ([T]), see for example 
[S], [H], and references therein. 

These results say that the large time asymptotic is given by the solution 
of a homogenized problem. It remains however to understand if this homog- 
enized problem leads to a non zero asymptotic averaged velocity. For this 
purpose we will consider the quantity 



which describes the asymptotic average displacement of the particle per unit 
time (provided the limit exists). 

Our first theorem states that the average asymptotic velocity is typically 
non zero. 

Theorem 1.1. The quantity I{b) is independent of the initial condition, and 
the set of b E with space average and time average equal to zero where 
I{h) 7^ is open and dense. 

Remark 1.1. By assuming that the space average (which may depend on 
time) and the time average (which may depend on space) are both zero we 
restrict the problem to a smaller set of possible drifts. One can prove a 
similar result with weaker constraints. Note also that it is well known that 
if b does not depend on time, then I{b) =0 (see for example [I3j). 

Remark 1.2. The theorem can be extended in various directions, for exam- 
ple by using different topologies on 6, by including a non-constant periodic 
diffusion coefficient, or by considering almost periodic time dependence. 

Another common model for molecular motor is the two states model which 
describes the time evolution of two non negative function pi and p2- In this 
model, the "molecule" can be in two states: 1 or 2 which have different 
interaction with the landscape described by the drift. We denote by pi(t,x) 



= lim - 

t— >oo t 




(2) 
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the probability to find the molecule at site x ate time t in state 1, and 
similarly for p2- We refer to [13] for more details. The evolution equations 
are given by 



dtpi = djDdxPi + hi{x)pi) -U1P1 + Z/2P2 



dtp2 = d^(Dd^p2 + h2{x)p2) + vipi - Z/2P2 



(3) 



where D, vi and z/2 are positive constants, hi and 62 are periodic functions 
of X of period L the last two with average zero. 

The asymptotic average displacement per unit time of the particle is now 
defined by 

/(z/i, z/2, 61, &2) = hm \ I X (pi(t,x) + p2{t,x)) dx . 

t^oo t J 

We have the equivalent of Theorem ll.il As before, we assume that |x| (P1+P2) 
is integrable. 

Theorem 1.2. For any constants ui > and z/2 > 0, /(z/i, z/2, &i, 62) is 
independent of the initial condition, and the set of hi and 62 G with space 
average equal to zero where /(z/i, z/2, &i, 62) 7^ open and dense. 

Another model of a particle interacting with a thermal bath is given by 
the Langevin equation 

dx = V dt , ■ 

dv = (-7 V + F{t, x) /m) dt + a dWt ^ ^ 

where m is the mass of the particle, 7 > the friction coefficient, F{t,x) 
the force, Wt the Brownian motion and a = V2D where D is the diffusion 
coefficient. We refer to [6] and [lO] for more details. 

For the time evolution of the probability density f{t, v, x) of the position 
and velocity of the particle one gets the so called Kramers equation 

dtf = -V dj + [{-IV - F{t, x)/m)f] + ^dlf . (5) 

We refer to [6] and references therein for more details on these equations. 
By changing scales, we can assume that m = 1 and D = 1 and we will only 
consider this situation below. Moreover we will assume as before that F{t, x) 
is periodic of period T in time, L in space and with zero average in space 
and time. We can now define the average asymptotic displacement per unit 
time by 

1(7, F) = lim - / dr /(r, f , x) dv dx . (6) 

t — ^00 t 
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As for the previous models the average asymptotic velocity is typically 
non zero. As usual, we denote by H^{dvdx) the Sobolev space of square 
integrable functions of x and v with square integrable gradient. 

Theorem 1.3. For 7 > 0, 1(7, F) is independent of the initial condition, 
the set of F & with space average and time average equal to zero where 
7(7, F) ^ is open and dense. 

One can also consider a situation where the particle can be in two states 
which interact differently with the landscape. This leads to the following 
system of Kramers equation. 



dtfi = ^d^fi - vdji + 
dtf2 = \dlf2 - vd^f2 + dy 



{^v-F,{x))f, 
(7^-F2(a;))/2, 



+ i^i/i - 1^2/2 ■ 



(7) 



In this equation, Fi and F2 are two periodic functions representing the 
different interaction forces between the two states of the particle and the 
substrate. The positive constants i/i and 1/2 are the transition rates between 
the two states. The non negative functions fi and /2 are the probability 
densities of being in state one and two respectively. The total probability 
density of the particle is the function fi + /2 which is normahsed to one. The 
asymptotic displacement per unit time for this model is given by 



7(7, Fi, F2, ui, V2) = hm - j ds j j V (/i(s, v, x) + /2(s, v, x)) dv dx , (8) 
and we will prove the following result 

Theorem 1.4. For 7 > 0, z/i > and 1/2 > 0, Fi, F2, ui, ^2) is inde- 
pendent of the initial condition, and the set of Fi and F2 € with space 
average equal to zero where 1(7, Fi, F2, i^i, 1^2) ^ is open and dense. 

2 Elimination of the spatial average. 

Before we start with the proof of the Theorems, we first show that the result 
does not depend on the spatial average of the drift b. 

Proposition 2.1. Assume b has space time average zero, namely 
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Then the drift b given by 

1 

b(t,x) = b(t,x + a{t)) — — / b{t,y)dy 

L Jo 

where ^ ^ 

a{t) = ds b{s,y)dy . 

is periodic of period T in time and of period L in x. This drift has zero space 
average and leads to the same asymptotic displacement per unit time. 

Proof. Note first that since the space time average of b is zero, the function 
a is periodic of period T. Let w be a solution of ([T]), and define the function 

v{t, x) = u{t, X + a{t)) . 

An easy computation leads to 

dtv = (^(^xV + b{t, x)v 

Since a{t) is periodic and bounded we have by a simple change of variable 

If If 

lim - / xu(t,x)dx = lim - / xv(t,x)dx . 

t— ►oo t J t^oo t J 

□ 



3 Proof of Theorems 11.1 



We start by giving a more convenient expression for the asymptotic velocity 
I{b). Using ([2]) and ([I]), and integrating by parts we get 

I{b) = lim - I ds I xdsu{s,x)dx = — lim - [ ds [ b{s,x)u{s,x)dx . 

t^OO t Jq J t^OO t Jq J 

Since b is periodic in x, of period L, we can write 

b{s,x) u{s,x)dx = / b{s,x) Uper{s,x)dx 



where 



^per 



[s, x) = u{s, X + nL) 



is a periodic function of x of period L. Note that since b is periodic of period 
L, Mper satisfies also equation ([T]). We now have 



1 /■* 

I{h) = — hm - ds b{s, x) Uper{s, x)dx . (9) 
Jo Jo 

Since the system is non autonomous, although periodic in time, we can only 
expect that when t tends to infinity, the function Uper{t, x) tends to a periodic 
function Wb(t,x) of t and x. Let Wb be the solution of equation ([T]) periodic 
in space and time and with an integral (over [0,L]) equal to one. It can 
be expected (see [13], [3], [TT]) that the asymptotic average displacement is 
given by 

/(6) = -1 / / b{t,x)wb{t,x)dtdx . (10) 
T Jo Jo 

In order to give a rigorous proof of existence of the function Wbit, x) and 
of the above relation, we introduce a new time. We consider the operator C 
given by 

Cw = —dsW + dx ^-dxW + b(s, x)w 

acting in a suitable domain dense in the space L^^^{dsdx) of integrable 
functions which are periodic in s and x of periods T and L respectively . 
This operator is the generator of the diffusion on the two dimensional torus 
([0,T] X [0, L] with the suitable identifications) associated to the stochastic 
differential equation 

ds — dt / 1 1 \ 

dx = -b{s,x)dt + dWt ^ ' 

where Wt is the standard Brownian motion (see [H]). We can now establish 
the following result. 

Proposition 3.1. The diffusion l[Tl\) has a unique invariant probability mea- 
sure with density Wb{s,x). This function is strictly positive. It is periodic of 
period T in s and of period L in x and satisfies equation (Qp, and it is the 
only such solution. The semi group with generator C associated to the dif- 
fusion ( fIT]) is compact and strongly continuous. The peripheral spectrum of 
its generator is composed of the simple eigenvalue zero (with eigenvector Wb). 
In particular, for any function v G L^^^{ds dx) , we have in the topology of 
Ij^^^{ds dx) 

lim e'^^f = Wb I I v(s, x) ds dx . 
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This kind of results is well known, we refer to for an exposition and 
further references. We can now establish the relation between (Q, and ffTUl) . 

Proposition 3.2. Let vq > be a periodic function of period L in x of 
integral one. Denote by v{t,x) the solution of ^\) which is periodic of period 
L in X with initial condition Vq. Then 

lim - I ds j b{s, x) v{s, x)dx = / / b{t,x)wb{t,x)dtdx . 



t^oo t In .In T 



^0 



Proof. In order to apply Proposition (13. ip . we consider the operator given 
by 

-Co = dx I -dxU + x)u 



2 

and observe that if w G L^^^^dsdx), we have 

(e^^w) (s, x) = {e^^°w{s - r, ■ )) (x) , (12) 
and in particular for any integer n, we get 

(e"^^«;) (s,x) = (e"^^«u;(s, ■ )) (x) , 

since w is of period T in s. 

We now take for w the function w{s, x) = v{s, x) for < s < T. Although 
this w may have a jump at s = T, we can consider it as a function in 
Lp^j.{dsdx). We observe that if W{t, s,x) is a solution of 

d^W = -dsW + dx (IdxW + b{s, x)W 



2 

then for each fixed sq, the function hs^^r, x) = W{t, sq + r, x) is a solution of 
drhso = dx (^^dxhso + b{so + r, x)/^^,, 

Therefore, by the uniqueness of the solution of ([T]), we have for any t > 
(taking so = t-T[t/T]) 

v{t,x) = (e[*/^l^^w7) (t - [t/L]T,x) . 

The proposition follows by applying Proposition 13.11 

□ 

The following proposition is the other main step in the proof of Theorem 



7 



Proposition 3.3. The function b ^ I{b) is (real) analytic in the Banach 
space . 

By this we mean (see [I2]) tliat the function is C°°, and around any point 
b E there is a small ball where the Taylor series converges to the function. 

Proof. We will establish that the map b \—>- Wb is real analytic in Lp^j.{ds dx) . 
For this purpose, we first establish that the operator A defined by 

Av = dx{bv) = bdxV + d^bv 

is relatively bounded with respect to 

and with relative bound zero (see [9] for the definition). This is obvious for 
the operator of multiplication by d^b which is bounded, and since b is bounded 
it is enough to derive the result for the operator d^- We will show that there 



I L^^^{ds dx) 



is a constant C > such that for any A > 0, [|9a.-RA| 

where Rx is the resolvent of C. In other words, we will show that for any 
A > 

< — 

Analogously to formula ( |T2l) we have for any w G L^^j,(^ds dx) 



a. / e-^^e^^dr 



(13) 



dx 



rC' 



W ) (s, x) 



dxgT{x,y)w{s - T,y)dy 



where gr{x,y) is the heat kernel on the circle of length L. We now observe 
that if n is an integer with > 2, we have (since x G [0, L]) 



sup 

ye[0,L] Jo 



\x — y — riL 



r3/2 

From the explicit expression 



.^x-y-nLf/{2r)^^ < ^(1)- 



-n2/(4T) 



gr{x. 



^-{x-y-nLf/{lT) 



it follows easily that 



sup 

j/G[0,L] Jo 



dxgr{x,y) dx sup / 

, l<^ye[o,L] Jo 



\x — y — nL I 



-3/2^2 



-{x-y-nLf/(2T)^ 



X 



|n|>2 ^ 



Therefore, we get 

9. fe^^'w; 



^ __\W\\Ll,,{dsdx) 



Lj,^j.{ds dx) 

Multiplying by e~^'^ and integrating over r we get the estimate (fT3ll which 
implies immediately that A is relatively bounded with respect to C with 
relative bound zero. Since the eigenvalue of £ is simple and isolated, 
the proposition follows from analytic perturbation theory for holomorphic 
families of type (A) (see [9J). □ 

In order to prove Theorem ll.il we now establish that / is non trivial near 
the origin. 

Lemma 3.4. DIq = and 

pqL^T^ . ,2 



P,<1 



•2 



Proof. The first statement is immediate from formula (ITOl) . For 6 = 0, we 
have obviously Wq = 1/L. For b small we use Taylor expansion. We have for 
b = eB 

As we just explained, w° = 1/L and we have for the equation 

Moreover, must have space time average zero. This equation can be solved 
using Fourier series in time and space. Namely if 



we obtain the equation 

2mp ^ ^ 2-niq 

or in other words for any {p, q) ^ {0, 0) 

9 



Note in particular that the denominator does not vanish except for p = g = 0. 
Using the Plancherel formula we can now estimate /(e-B). We have 



Therefore 

Since B is real, we have Bp^g = B^p^^q and this can also be written 

PA 

This finishes the proof of the Lemma. □ 



To prove Theorem ll.H we observe that since / is continuous, the subset 
of where it does not vanish is open. If this set is not dense, the zero 
set of / contains a ball. However since I is real analytic and is pathwise 
connected we conclude that in that case / should vanish identically contra- 
dicting Lemma 13.41 We refer to [12] for more properties of the zero set of 
analytic functions in Banach spaces. 

Remark 3.1. We observe that D^I is a non definite quadratic form. This 
leaves the possibility of having non zero drifts h (with space and time average 
equal to zero) satisfying I{h) = 0. Let bi and 62 be such that > and 
7(62) < 0. Such bi and 62 exist, one can for example take them of small 
enough norm and use Lemma 13. 4[ Moreover one can assume that 62 ^ R&i • 
Otherwise, by the continuity of J, one can perturb slightly 62 such that this 
relation does not hold anymore, but /(&2) is still negative. One now considers 
the function if (a) = J ((1—0)61+062) • This function is continuous, it satisfies 
ip{l) > and ip{0) < 0, hence it should vanish at least at one point ao g]0, 1[. 
At this point we have bo = {1 — q;o)6i + 00^2 7^ and /(60) = 0, and bo is a 
non trivial periodic function with vanishing space average and time average. 



4 Proof of Theorem II. 2L 

In this section we discuss the model with two components ([3]). As before, we 
arrive at the formula 

J(z/i,z/2,6i,62) = - lim J / / (bi{s,x) pi{s,x) + b2{s,x) p2{s,x)) dsdx , 

t Jo Jo 
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where pi and p2 are solutions of periodic in x of period L. We denote 
by Lpgj.((ia;) the space of integrable periodic functions of period L in x with 
value in R^. The norm is the sum of the norms of the components. 

Proposition 4.1. The semi-group defined by ^ is compact in Lp^j.{dx). It 
is positivity preserving, and its peripheral spectrum is the simple eigenvalue 
one. The corresponding eigenvector can be chosen positive with dense support 
and normalised, it depends analytically on bi, and 62- 

Proof. We introduce the three generators 

\ P2 ) V (^x{DdxP2 + b2{x)p2) + VlPl - h'2P2 ) ' 



\P2 J V Ddlp2 J \P2 J \ d,{Dd^p2 + b2{x)p2) J 

This operator Cq is the infinitesimal generator of a strongly continuous 
bounded and compact semi-group in Lp^^^dx). It is easy to verify that the 
operator A = C — Cq is Cq relatively compact and £0 relatively bounded with 
relative bound zero (see |9j). Therefore C is also the infinitesimal generator 
a strongly continuous and compact semi-group in Lp^j.{dx), and similarly for 

A. 

The semi-group e*^ is positivity improving (see [2]). Indeed, let M be the 
matrix 

-Ul V2 



M 

yi -V2 

It is easy to verify for example by direct computation, that the matrix e*^ 
has strictly positive entries for any t > 0. Moreover, for any t > 0, we have 
gtAf ^ e~**^'^^"^'^^''ld, where the inequality holds for each entry. It immediately 
follows from the Trotter product formula (see [2] ) that for each x and y in 
the circle, and any t > 0, we have 



{x,y) > e-'^'''^'''^e'^'ix,y) 



again in the sense that the inequality holds between all the entries. Since for 
each X, y and t > the diagonal elements of e*^^ (x, y) are strictly positive 
(see for example [15j), we conclude that the matrix valued kernel e*^(x, ?/) 
has non negative entries and strictly positive entries on the diagonal. 

Since the sum of the integrals of the two components of an element of 
Lpgj.{dx) is preserved by the semi-group e*"^, it follows that this semi group 
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has norm one in L^^^{dx). It then follows by classical arguments that is a 
simple isolated eigenvalue of the generator and there is no other eigenvalue 
with vanishing real part. 

The analyticity follows from the uniqueness and simplicity of the eigen- 
value as in the proof of proposition (13. 3p . 

□ 

We denote by wi and W2 the two (non negative) components of the sta- 
tionary solution of the system ([3]) which are periodic of period L and nor- 
malised by 

[wi{x) + W2{x)~\ dx = 1 . 

Note that wi and W2 depend on the constants z/i, z/2, and the functions bi, 
and 62- 

It follows immediately from Proposition 14.11 that the average asymptotic 
velocity is given by 

/(z/i, z/2, 61, 62) = / [bi{x)wi{x) + b2{x)w2ix)]dx . (14) 
Jo 

Since the function /(z^i, z/2, ^2) is analytic in (bi, 62), to prove that it is 
non trivial, we look at the successive differentials at the origin. 

Proposition 4.2. DIq = and for any (61,62), -D^/o((6i, 62), (61, 62)) = 0. 

Proof. The first result is trivial. For the second result one uses perturbation 
theory as before in the Fourier decomposition. We get with a = Air'^D/L'^ 

—an^ — vi V2 \ _ 2% in fi'2bi{n] 

z/i -an'^-U2j \wl{n)J ~ L{ui + U2) \'^ib2{n] 

Some easy computations using Plancherel identity lead to 

D^/o((6i, 62), (61, 62)) =-y: ^ 2TZ , — Y 

X {an^ + z/i)z/26i(n)6i(n) + [an^ + i'2)vj)2{n)b2{n) 

-Fz/i z/2 (62(^)61(72) +61(^)62(72)) 
since 61 and 62 are real (61(72) = bi{—n)). □ 
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This result suggests to look at the third differential at the origin which 
turns out to be a rather involved cubic expression. In order to show that 
the function / is non trivial, it is enough to find a particular pair (&i,62) 
such that D^Io{{bi, 62), {bi, 62), (&i, ^2)) 7^ 0. This was done using a symbolic 
manipulation program (Maxima). We found that for L = 271, D = 1, bi{x) = 
cos(2x) and &2(a^) = cos(x), one gets 

D%{{b,M),ibM,ib,M))- ^M-2-2u, + l) 



4(z/2 + Ui){u2 + Ui + l)2(z/2 + Z/i + 4) ' 

Theorem 11.21 follows as before. 



4.1 Proof of Theorem 11.3 

As in the previous section, we can introduce the periodised function (in x) 

fit, v,x) = f{t, v,x + nL) . 

n 

This function is periodic of period L and satisfies also equation ([5]). We get 
also li^-^^F) by replacing / by f(t,v,x) in equation ([6]) and integrating only 
on one period. From now on we will work with this periodised function and 
denote it by / by abuse of notation. 

We now introduce a stochastic differential equation on [0, T] x [0, L] x R 
with periodic boundary conditions in the first two variables s and x. This 
differential equation is given by 

ds = dt 

dx = vdt (15) 
dv = —yvdt + F{s, x)dt + dWt ■ 

To this diffusion is associated the infinitesimal generator C given by 

Cw = —dgW — V dxW + [(7f — F{t, x))'w\ + i^dlw . 

We denote by the space L'^{e^'"^ds dv dx) of functions periodic in s of 
period T and periodic in x of period L. Using an space instead of an L} 
space is useful in proving analyticity. 

We can now establish the following result. 

Proposition 4.3. The diffusion semi-group defined by 173]) in SS is compact 
and the kernel has dense support. It is mixing and has a unique invariant 
probability measure (in SS) with density f{s,v,x). This function is strictly 
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positive, satisfies equation and it is the only such solution. In particular, 
for any function w E we have in the topology of ^ 



Mm e'^^w = f / / / w{s,v, x) ds dv dx . 
Jo Jo J 

The function f is real analytic in F E . 

Proof. Instead of working in the space ^ , we can work in the space LF'{ds dv dx) 
by using the isomorphism given by the multiphcation by function e"^" In 
that space we obtain the new generator Cp given by 

C'pg = -dsg -vd^g- Fd^g + F-fvg —g + -g + -d^g . 

Using integration by parts, it is easy to verify that 

^ g {C'pg^ds dvdx = 



+ 7^^ Ifi-P - Y^^kP - ^ j \dv9f^ dsdvdx 

< (^■^ + 11-^11^) J \g\'^ ds dv dx — J (^v'^\g\'^ + '^\dvg\'^^ dsdvdx. (16) 

We see immediately that —Cp is quasi accretive (see [S] for the definition 
and properties of the semi-groups generated by these operators). 
Let gt = e^^fgo, using several integrations by parts one gets easily 



/y^ /'||2 ||||2 f \ \2 

< — — / "W^ (isc/fc/a; + (7 + / \gt\ dsdvdx. 

We obtain immediately 

and 

tdr [ v' \gr\' dsdvdx < e(^+^^'ll^l'-)* H^oH' ■ (18) 
Jo J 7 
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Similarly, we get 

9t I kiT ds dv dx 



4 J Fv\gt\ dsdvdx + 2'j J Fv \gt\ dsdvdx — 'j J v \gt\ dsdvdx 
+7 / v"^ \gtf ds dv dx + 2 / IgA'^ ds dv dx — v^ Id^gtl'^ ds dv dx 



< C(7, ||F||oo) J \gt\ dsdvdx (19) 

where C{-y, ||-F||oo) is a constant independent of gt- For t > fixed, we deduce 
from ffTSj) that there exists ^(t) e [0,t[ such that 

/ v^ \g^^t)\'^ ds dv dx < e(^+^^'ll^ll-)* \\go\\'^ ■ 

Using (14.11) and (|T71) we get (for any t > 0) 

J v^ \gt\'^ ds dv dx = J v'^ \g^[t)\'^ ds dv dx + J dr Or J v"^ \gr\'^ ds dv dx 



'm 

|2 



1^0 I, 

In other words, for any t > the image of the unit ball by the semi-group 
is equi-integrable at infinity in v. Compactness follows immediately using 
hypoelliptic estimates (see [IS])- 

From the standard control arguments (see [I^) applied to the diffusion 
ffT^ . we obtain that the kernel of the semi-group has dense support. We now 

2 

observe that integration of a function against e~'^'" is preserved by the semi- 
group evolution. This implies by standard arguments that the spectral radius 
of the semi-group is one, the invariant density is unique, and exponential 
mixing holds (see [13]). Finally, as in equation fl4.ll) . for g in the domain of 
C we have 



^ J 9 {^'f 9) dv dx 

> — 2||F||^j J \g\'^ ds dv dx — J (^j^v'^\g\'^ + -\dvg\'^^ dsdvdx 
This implies for any A > 

j \d^g --ivg\^ dsdvdx < '^\\C'g\\l+ (4 A + 3 7 + 4 ||f||^) ||^||2 . 

15 



In other words, the operator F (d^ — 7f ) is relatively bounded with respect 
to £q with relative bound zero and this implies the analyticity (see |9]). □ 

We can now complete the proof of Theorem ll.3[ 

Proof, (of Theorem 11.31) . 

Repeating the argument in the proof of Proposition [3]2l using Proposition 
14.31 we get 

lim - / / / V f(T,v,x)dT dv dx = — [ [ [ F(t,x) f(T,v,x)dT dv dx , 

We can now use perturbation theory to compute the right hand side near 
F = 0. For this purpose, it is convenient to fix a function G periodic in 
space and time and with zero average, and consider F = eC with e small. 
Since / is analytic by Proposition 14.31 we can write 

/ = /o + e/i + e72 + ^(l)e' 

where 

and for n > 1 the /„ are functions of integral zero, periodic in time of period 
T, defined recursively by 

dtfn + Vdjn - dy{'yvfn) - ^dlfn = -G{t, x)d^ fn-1 ■ 

We get immediately 

mG(r, x) /o(t, V, x)dT dvdx = , 

since /o is independent of x and t and G has average zero. Therefore we now 
have to look at the next order perturbation, namely the second order in e for 
the average velocity. In other words, we have 



t— too t 
2 pT pL 




lim - / / / V f{T^v, x)dT dv dx 



'OJO 



e 

jT Jo 
We first have to solve 




G{r,x) fi{T,v,x)dTdvdx + ^(l)e^ (20) 



dtfi + vdji - d^{-fvfi) - ^dlfi = -G{t,x)dJo 
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to get /i. For this purpose, we use Fourier transform in all variables (recall 
that t and x are periodic variables). We will denote by /i,p,g(^) the Fourier 
transform of /i (and similarly for other functions), namely 



pt/T 2TTiqx/L I ikv 



P,<1 



e*'Vi,p,,(A;) dk. 



We get 



2Tcip k"^ 



fl,p,q + ik 



27rg\ d 



L J dk 



ikG. 



271 L 



hi £-^=7(47) 



(21) 



We now observe that equation (14. ip . only involves the integral of /i with 
respect to v since G does not depend on v. Therefore we need only to 
compute /i,p,q(0). Let hp^g{k) be the function 



1 - 



'jkL 



27r^g-^/{7-^L'^)+27rip/{7T)-l 



For — 27r|g|/(7L) < k < 2'K\q\/{'jL), this function is a solution of 

27Tq 



2TTip , k'^\ d 

i^p,q\f^) 



T 



dk 



^k 



hp,q{k^ 



. 



For g > 0, multiplying ( |2T]) by hp^q{k) and integrating over k from to 
27rg/ (jL), we get 

/l,p,ij(0) CrpqYpq 

where 



27rg/(7L) 



,TTqk/{-y'^L) 



1 - 



'-fkL 



2nq 



2iT^q^/{-y^L'^)+2TTip/{-yT)-l 



kdk . 



Note that since q ^ 0, the integral is convergent. For g < 0, one gets a 
similar result, namely 



^nqk/("f^L) 



27rg/(7L) 



'ykL 



2Txq 



2-K^q^ /{r^L^)+2-Kip/{jT)-l 



k dk , 



and it is easy to verify that Tp^q = r_p^_q. We now have 




'0 ^0 



G(r, x) /i(r, V, x)dT dv dx = Gp^q Tp^q Gp^, 



p,<i 
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Since G{t, x) is real, we have Gp^q = G-p-g, and therefore 
G'(r, x) /i(r, t>,x)c?r(it>(ix = [Gp ' 

We observe that for g > 0, 




'0 ^0 



r + r 1 



27rg/(7L) 



27r2g2/(73L2)-l 



X sin ((27rp/(7T)) log(l - -fkL/{2TTq))) kdk , 
= e2-''?'«/(^'^')(l - «)2-^'?V(73i-^)-i sin ((27rp/(7T)) log(l - u)) u du . 

Note that this quantity is equal to zero for p = 0, and it is odd in p. It can 
be expressed in terms of degenerate hypergeometric functions. We now have 
to prove that for p 7^ 0, this quantity is not zero, at least for one pair of 
integers {p,q)- For this purpose, we will consider q large. To alleviate the 
notation, we consider the asymptotic behavior for a > large and /3 G R 
fixed of the integral 







Using steepest descent at the critical point u = (see 0), one gets 



We apply this result with a = 27r^g^/ (7^^^) and /3 = 27ip/ (7T), and conclude 
that for q large enough, [Tpg + T_p _gj /2 7^ as required. □ 

The proof of Theorem 11.31 is finished as before. 



4.2 Proof of Theorem 11.4 

The scheme of the proof is similar to the proofs of the previous Theorems. 
We only sketch the argument except for some particular points. We assume 
7 > 0, z/i > and 1^2 > 0. One starts by reducing the problem to a periodic 
boundary condition in x. The key result is the analog of Proposition 14. 1[ 
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Proposition 4.4. The semi-group defined by ^ is compact in L'^ie^'" dvdx) 
(the X variable being on the circle of length L). It is positivity preserving 
positivity improving on the diagonal. Its peripheral spectrum is the simple 
eigenvalue one. The corresponding eigenvector (/i, can be chosen positive 
with dense domain and normalised (fi + /2 of integral one). The functions 
fi and /2 depend analytically on Fi, and F2. 



Proof. The proof is very similar to tlie proof of Proposition 14.11 and we only 
sketch the details for the different points. 

□ 

It follows as before from the evolution equation that 

/(7, Fi, F2, 1^1, V2) = j {fiFi + /2F2) dv dx . 

For fixed 7 > 0, z^i > and z/2 > 0, this quantity is real analytic in (Fi, F2), 
and to check that it is non trivial we investigate its behaviour near the 
origin. For this purpose, we set (^1,^2) = e{Gi,G2) with Gi and G2 two 
functions, periodic of period L and with zero average. We now develop 
7(7, eGi, eG2, 1^1, 1^2) in series of e. As for model ©, the terms of order e and 
of order vanish. One then has to find a pair (Gi, G2) such that the term 
of order does not vanish. The computations are rather tedious and will be 
detailed in the appendix. 



A Proof of Theorem 11.4 



We define the vectors 



/=( f^] and F ( 



We then have 



7(7, Fi, F2, z/i, z/2) = jj F*f dvdx. 
The equation for the stationary solution of ([7]) can be written 

]-dlf- vdj+ [jvf~] + Mf = B{x) dj 



[22) 



where as before 



M 



-I'l V2 
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and 

Fi(x) 



B x) 



Foix 



Without loss of generality, we can rescale v and assume 7 = 1, replacing z/i 
and U2 by ui/'y and 1^2/7 respectively, L by L7~^/^, Fi and F2 by F\j ^ and 
-^2 / a/7 respectively. 

For reasons that will become obvious later on, it is more convenient to 
consider instead of / the vector 

Let .if be the operator defined by 

££h[v) = dlh{v) - v^h{v) + h , 
then the equation ([22]) for the stationary solution can be written 

- 2vd^tp + 2Mtp = 2B{x) [d^^ - vip) , (23) 

and we have 

j(7,Fi,F2,z/i,z/2) = jj F.^fe-^"'^ dvdx, (24) 
with normalised by 

{1P1+1P2) e"'''/^ dvdx = l. (25) 



For F = eG we have the expansion 

^=^o + e7/^^ + eV^ + ^(e^) (26) 

where 



and the vectors ip ^ are obtained recursively by solving 

^iP^^^ - 2vd^^^^^ + 2MV^^+^ = 2C{x) {d^ip^ - vip^) , 

where 

_ ( G,{x) 
^^^^ ~ V G2{x) 
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Note that ^ satisfies 

^^p^ - 2vd^i^^ + 2M^^ = 



and the normahsation condition (l25l) . while for j > 1, should satisfy the 
normalisation condition 



(iji + iji) e-" /2 dvdx = 0. (27) 
Corresponding to the expansion fl2U|) . we have the expansion 

3 

/(l,eGi,eG2,z/i,z/2) =$^e"/n(G'i,G2,z/i,i/2) + ^(e") , 



n=l 

where 



We have immediately Ji (Gi, G2, z^i, 2^2) = since the averages of Gi and G2 
over X vanish and ip ° does not depend on x. In order to compute Ji, we need 
to compute ip ^ which solves 

^if^ - 2vd^^^ + 2M^^ = ve-'''/^H{x) 

where 



H(x) 



z/2Gi(x) 



Since if is periodic of period L, we can decompose this equation in Fourier 
series and get 

We now observe that 

^ _ Ampv ^ 2M = e^-^pS"/^ (^^ - 1^ + 2M^ g-^-^f ^"/^ . 

This implies 

(/I 2 2 \ ~i 

^_Z^ + 2Mj e-^^"'^-'^ (ve-^'/^Hip)^ . (28) 
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The spectrum of the self adjoint operator in L'^idv) is well known (see for 
example [T] or [3]). The eigenvalues are the numbers — 2n with n G N*, and 
the corresponding normalised eigenvectors are the functions 



en[v 



( ^ \n 9— n/2 Jn 

Hjv) e-''/2 = ^-^^ — ^-'^^ 



7rV4(n!)V2 ^ 7^1/4(^1)1/2 dv 

where the Hn are the Hermite polynomials. We have 

271 ip 



g -2nipd^/L (^Q-v'/2 



L 



and we decompose this function on the Hermite basis (e„). It follows easily 
using integration by parts that 



en{v) V 



277 



-t. V2g27r j p g2 7r2 p2 /i2 



n 2 7r2p2//,2 



-1)'^ e 



2-n/2 



7ri/4 [n\y/2 



,2 7r2p2/L2_ 



27Tip 



Jn 

^2^ipv/L ^g-1.2^^ 



2-n/2 



2-2/L2 



271 ip 
L 

,2-"/27,l/4 



vri/4 (n!)i/2 
'27cip 



X 



+ n 



L 



-v^^2nipv/L 



(n!)i/2 



(n + 27rV/L2) 



271 ip 



n — 1 



This formula holds for any p 7^ and any n > 0. For p = 0, the formula 
holds for any n > 1 (with the convention O" = 1). Finally, for n = and 
p = the result is zero. We now have after some simple linear algebra 



^ - + 2M 



-1 



A7['^p^ 

~1X 



enH{p)) (v) = en{v) { -2n 
en{v) 



-1 



+ 2M H{p) 



2{n + 27rV/^^) [n + vi + V2 + 27i^p'^ / L"^) 
-V2 —n — 27\'^p^ jL^ —V2 



H{p) 



-ui-n- 27rV/L2 
We can write using Fourier series 
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p n ^ ^ 

2(n + i/i + z/2 + 27rV/^^) ("-O^^^ 



Since the operator id^ is self adjoint, we get 



By a computation similar to the above one, we get 

After some simple algebra, we obtain 

/2(Gi,G2, 21/1, 2^/2) 



P 



where 



2 V- 2.^pVi^2-^ /27r^py"-^ 

+ n! V L ; (n + + i/2 + 2vrV/^^) 



= Gi(p)G'i(p)z/2(z/2 + n + 27rV/L2) +G2(p)G'2(p)z/i(i/i + n + 27rV/L2) 



+ Z/lZ/2(Gi(p)G2(j9)+G2(j9)Gi(p)) . 



Since Gi and G2 are real, we have Gi{p) = Gi{—p) and G2{p) = G2{—p)- 
This implies immediately that R{p) is even in p and therefore /2 (Gi, G2, z/i, z/2) 
0. A similar computation (involving only p = 0) shows that the normalisa- 
tion condition (p71) holds. We now need to compute /3(Gi, G2, z^i, 2^2) , and 
we recall that it is given by 

h{G^,G2,u^,U2)=J2'^)^ [ ^\v,p) e'^" dv (29) 

where 
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and where J{v,p) is the Fourier series of the vector 

J{v, x) = 2C(x) (dyif^v, x) - vif^v, x)^ . 
As above, using that the operator idy is self adjoint, we get 



p 

We can now compute everything in terms of series of Hermite coefficients. 
Denoting by C{p) the Fourier coefficients of the matrix C{x) we get 



Q 

Q 

We have already computed (for n > if p 7^ and for p = if n > 0) 



^ (n!)V2 L ; 2(n + i/i + i/2 + 27rVA') 



—i'2 — n — 2n'^q^/L'^ —i>2 

—vi —vi — n — 27r^g^/L 

By a similar computation left to the reader, we get 



2.2 /r2 \H{q) (30) 
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e 

m 



(m!)i/2 



X ■ 



2(m + 27r2p2/^^) (m + z/i + z/2 + 
— z/2 — m — 27r2p2/^2 

It finally remains to compute 



(31) 



er,^(w)e-'"*(P-'?)^"/^(9, -v + 2mq/L)en{v)dv 
-V2y/n + 1 / e„(f)e„+i(t; - 2 7ri (p - (it; 



since 

(c^ ^ -V2Vn + 1 e„+i(t;) . 

We now compute the quantity 

ln,m{r) = J em{v)en{v -271 ir/L) dv 

2-(n+m)/2g2 tt^ r^/L^ 



a/tF (n! m!)^/2 



^/7^ (n! m!)^/2 

2-(n+m)/2g7r2rVi2 



-^/tt (n! m!)^/2 
When m < n, we obtain 



y g-(.-.^r/L) Hm{v) Hn{v-2 7rtr/L) dv 

J e~^^ Hm{v + nir/Lj Hn{v — vrir/L) (if . 



7n,„.(r) = ^^-^ \ ( - ^ I 



n!)V2 \ L J \ 

where denotes a Laguerre polynomial (see [4J, formula 7.377). Similarly, 
when m > n we get 

^ ^ _ (^!)l/22(>n-n)/2g.^rVL^ / ^ TT r ^^ ^ / 2rr^ 

1n,mV) — / iM/9 1 r ) 



(m!)V2 V ^ / "V 
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We can now use this result together with equations (jX]) and (jX]) into the 
expression (1^ . We get 

h{Gi, 02,1^1,1^2) = 

--rr~, — V ( — T^7n,m(p-g) - V^iTT7n+i,m(p-g) ) 

L{l^l + l^2) V / 

V / n,m,p,q ^ ' 



(n!)i/2 V ^ (m!)i/2 ^ l 



(m + {m + Ui + U2 + 2iT^p'^/L^) {n + Ui + U2 + 27r^q^/L^) 

where 

c/ X ( ( V2 + m + 27r V/^^ '^i \ 

^^^'^^=U ^2 ^1 + m + 27r V/^^ J ^^^T 

Gi{p-q) \ / 1^2 + n + 27r2gVL2 z/2 

G2{p-q) J \ ui z/i + n + 27r V/^' 



X 



^2Gi{q) 
T^iG2{q) 

A numerical simulation using for Gl and G2 real Fourier polynomials of 
order two and random coefficients give a nonzero result for h{Gi,G2,l,2) 
and L = 10. 
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